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An  investigation  of  stresses,  arising  from  a  thermal  field,  in 
electronic  packages  comprised  of  a  ceramic  chip  attached  to  a 
composite  board  by  a  solder  joint  was  undertaken.  Two  types  of 
solder  joint  connections,  a  leaded  device  and  an  unleaded  device, 
were  analyzed.  The  stresses  result  from  the  mismatch  of  material 
properties  of  the  tri-assembly  system.  Procedures  are  shown  for  the 
analyses  of  the  two  types  of  solder  connections.  The  analyses 
showed  that  the  thermoelastic  stresses  for  the  stiff  unleaded 
connection  are  significantly  greater,  by  a  factor  of  20  for  a  the 
particular  application  undertaken  in  this  report,  than  the 
thermoelastic   stresses   in   the  flexible   leaded   connection. 


Introduction 

Electronic  chip  packages  are  comprised  of  several  components  with 
different  material  properties  and  geometries.  One  of  the  most  basic 
electronic  packages  is  the  tri-assembly  system  consisting  of  a  chip 
attached  to  a  board  by  a  solder  connection.  During  the  fabrication 
of  this  simple  unit,  the  soldering  process  subjects  the  device  to 
a  thermal  field.  Since  the  three  components  of  the  device  have 
different  material  properties,  in  particular  their  Young's  modulii, 
Poissons  ratios,  and  thermal  coefficients  of  expansion,  the  mismatch 
in  these  properties  result  in  thermoelastic  stresses.  In  addition  to 
the  thermal  field  arising  from  fabrication,  thermal  fields  and  hence 
thermoelastic  stresses,  result  from  the  operation  of  the  device  itself. 

Not  only  is  the  mismatch  in  material  properties  important  in  the 
thermoelastic  stresses  that  develop,  but  equally  important  is  the 
manner  in  which  the  solder  connection  is  made.  For  a  given  set 
of  chip,  board  and  solder  properties,  an  important  question  'what 
is  the  best  way  to  make  the  solder  connection?'.  With  regard  to 
geometric  configuration,  there  are  two  main  types  of  connections, 
the  leaded  connection  and  the  unleaded  connection.  In  the  leaded 
connection,  the  chip  is  attached  to  the  board  by  a  thin  'lead'  as 
shown  in  Figure  1.  In  this  configuration,  the  chip  is  somewhat 
removed  from  the  board  its  attached  to.  In  contrast,  in  the 
unleaded  connection,  the  chip  is  attached  along  its  edges  to  the 
board  substrate  by  a  thin  layer  of  solder.  This  unleaded  joint 
arrangement  is  significantly  'stiffer'  than  the  leaded  connection,  and 
all  other  things  being  equal,  one  would  expect  that  the  stiffer 
unleaded  joint  would  result  in  larger  thermoelastic  stresses.  This 
investigation  seeks  to  determine  the  relative  magnitudes  of  the 
stresses  from  these  two  kinds  of  solder  joints.  It  should  be  noted 
however,  that  the  proximity  of  the  chip  to  the  board,  in  the 
unleaded  joint,  might  result  in  a  more  uniform  thermal  field  thereby 
reducing  the  thermal  gradients  and  thus  the  thermal  stresses.  The 
question  of  the  thermal  field  that  results  from  each  of  the  solder 
joints  is  not  addressed  in  this  investigation  which  only  considers 
the  effect  of  the  same  uniform  thermal  environment  on  each  of 
these    (leaded    and    unleaded)    solder  joints. 


Analysis    of    leaded    device 

Figure  1  shows  the  geometry  and  dimensions  of  the  leaded 
device  considered  in  this  study  which  consists  of  a  ceramic  chip, 
lead  wires  and  a  circuit  board.  The  cross-sectional  dimension  of  the 
lead  wire  is  given  in  Figure  2. 

Three  general  assumptions  are  made.  First,  the  whole  device  is 
under  a  thermal  cycling  load  of  +65°c  ~  -20°c  and  the  lead  wire  has 
an  initial  strain  from  soldering  process.  The  device  is  soldered  to 
the  circuit  board  at  a  temperature  of  +185c  and  this  tempature  is 
assumed  to  propagate  only  to  the  lead  wire  during  the  soldering. 
Another  assumption  is  on  the  deformation  of  the  device  that  the 
ceramic  chip  and  the  circuit  board  are  so  much  stiffer  than  the  leads 
so  all  the  deflections  due  to  the  thermal  loading  occur  only  to  the 
leads.  Finally,  the  device  is  assumed  to  be  geometrically  symmetric 
about  the  center  so  the  following  calculation  is  performed  only  for 
the  half  of  the  device. 

Material  properties  of  the  chip,  lead  and  the  glass-epoxy 
substrate  of  the  circuit  board  are  as  follows; 

Thermal  expansion  factor  of  the  chip,  cc_  =  6.5  10"6  ppm/=c 
Thermal  expansion  factor  of  the  lead,  aL  =  4.47  10-6  ppm/Dc 
Thermal  expansion  factor  of  the  substrate,  as  =  16.  10~6  ppm/°c 


Young's  modulus  of  the  lead  wire,  E  =  21.  106  psi 
Moments  of  inert 
=   1.77   10-10  in4 


Moments  of  inertia  of  the  lead  wire,      =  bh3/12   =  0.017  (0.005)3/12 


Calculation    of   thermal    mismatch 
At  +20cc; 


Initial    deformation    of    a    lead    wire    after    the    welding 
(from   +185°c  to  +20cc) 

A|N|T  =  a  AT(L)  =  4.47  1Q-6(-1  65c)(0.045)   =  -33.2   10"6  in 


At  +65°c; 


At   -20°c; 


Deformations  due  to  temperature   changes   from   +20°c  to 
+65°c 

Achip    =  acAT(Lc)  =  6-5  10-6(45°c)(0.26/2)  =  38.0   10-6  in 
Alead   =aLAT(LJ  =  AA1  "IO-6(45°c)(0.045)  =  9.1    10'6  in 
Asub    =asAT(Ls)  =  16.  10-6(45°c)(0.35/2)  =  126.0   10"6  in 

Mismatch  between  the  top(chip  and  lead)  and  the 
bottom    portion(substrate)    is 

A65°    =Atop-Abot 

"  <Achip  +  Alead  +  AINIt)  *  (Asut>) 

=   (38.0   10"6  +  9.1   10"6-  33.2   10"6)  -  (126.0  10'6) 

=   -112.1    10"6    in 


Deformations  due  to  temperature  changes  from  +20°c  to 
-20cc 

Achip  =  acAT(Lc)  =  6.5  10-6(-40cc)(0.26/2)  =  -33.8  10"6  in 
Alead  =  aLAT(LL)  =  4.47  10-6(-40°c)(0.045)  =  -8.05  10"6  in 
Asub    =ccsAT(Ls)  =  16-  10"6(-40cc)(0.35/2)   =  -112.    10"6  in 

Mismatch  between  the  top(chip  and  lead)  and  the 
bottom    portion(substrate)    is 

A-20*     =  Atop  '  Abot 

=  (Achip  +  Alead  +  AIN!T)  -  (Asub) 

=    (-33.8    10-6   -8.1    10'6-   33.2   10-6)   -   (-112.   10-6) 

=  36.9  10-6  in 


Since  A65c  is  much  greater  than  A_2q°  'n  its  magnitude,  it 

is   consider* 

calculation. 


s   considered   that  A65C   becomes   critical   for  the    stress 


Analogy   of  temperature   change   effect   and   support   settlement   effect 
It  can  be  seen  that  the  force  balance  from  the  temperature 


changes  are  the  same  as  that  from  the  support  settlement  shown   in 
Figure   3.    So   we    now   treat   the   displacement  A65o   from   the   above 

calculation  as  the  settlement  5. 

6  =  A65o 

To  find  the  force  P  cause  by  a  displacement  5  which  is  known 
from  the  previous  calculation,  we'll  break  the  frame  ABC  into  beam 
AB  and  beam  BC  as  in  Figure  4  and  5. 

Denoting  l_v  and  l_h  as  the  length  of  the  beam  AB  and  BC, 
respectively,  the  tip  displacements  and  the  tip  rotations  at  the  beam 
AB  due  to  an  applied  force  PB  and  a  moment  Mb  are  given   by  the 

formulas 

PLV3      MLV2 

5  =  ^eT"^eT  (1) 

_PLV2      ML^ 

eB-  2E|   -   El  (2) 

and  for  beam  BC  the  tip  rotations  can  be  expressed  by 

MLh 
0B=B  (3) 

Compatibility    requires   that   0B  at  the  above  equations  (2)  and 

(3)  be  the  same  so  by  equating  these  two  equations  an  expression  for 
M  can  be  found. 

PL2 
M-^^)  (4) 

Substituting    (4)    into    (1) 

PLV3      Lv  +  4Lh 
b~      El      12(LV  +  Lh) 

or  P  can  be  written  in  terms  of  5, 


SEI    12(LV  +  Lh) 
V  =  L  3       L     +  4L 


"V 


v  T  ^h 


_  112.1    10-6  (21.   1Q6)  (1.77  10-6)  12(0.028  +  0.045) 

(0.028)3  0.028  +  4(0.045)    "  °081      b 

Shear   stress   calculation 

The  shear  stresses  at  the  solder  joint  are  obtained  by  dividing  the  force 
by  footage   area   of  the   lead   wire   which    is   soldered   to   the   circuit 
board, 

_P  _  P  _      0.0811 
x"  A  "bl_"0. 017(0. 03)  "  159'°    PSI 


Analysis    of    unleaded    device 

Figure  6  shows  the  geometry  and  dimensions  of  the  unleaded 
device  which  consists  of  a  ceramic  chip  fixed  to  a  circuit  board 
by  a  thin  strip  of  solder.  In  the  present  application,  the  solder 
adhesive  layer  covers  only  a  small  area  along  the  edge  of  the 
chip. 

Two  models  of  the  problem  follows.  The  first  model,  developed 
by  Suhir  (1),  provides  an  analytical  determination  of  thermally- 
induced  stresses  in  a  tri-material  assembly.  The  second  model  is 
a  finite  element  model  which  utilizes  a  recently  developed  element 
which  maintains  inter-layer  continuity  along  the  chip-solder  and 
board-solder    interfaces. 

Suhir's    Model 

A  free  body  diagram  of  the  tri-material  assembly,  showing  the 
notation  used  hereafter,  is  shown  in  Figure  7.  The  Suhir  model 
is  applicable  to  assemblies  with  continuous  attachment  (a  =  I), 
as  well  as  non-continuous  attachment  (a  <  I).  The  development 
is  carried  through  by  invoking  equilibrium,  displacement 
compatibility,  and  thermoelastic  constitutive  relations.  The 
displacement  compatibility  relations  impose  the  condition  of  inter- 
layer  displacement  continuity.  For  a  detailed  development  of  the 
analytical  relations,  the  reader  is  referred  to  the  original  work. 
Here  a  brief  summary  of  the  development  is  presented,  followed 
by  a  sequential  list  of  the  calculations  that  lead  to  the  maximum 
interlayer    shearing    stress. 

Denoting    the    axial    displacements    at    the    top  and    bottom    of    the 

ith      layer      as      u+,(x)      =      u'^x)      respectively,  the      displacement 

compatibility  conditions  along  the  interfaces  of  layers  1-2  and 
layers    2-3    are: 

ul  (x)   =  ul  (x)  (5) 


u2'(x)    =  uUx)  (6) 

Axial  displacements  u'^x),  u+2(x),  u'2(x),  and  u+3(x)  are  obtained 
by  adding  the  individual  contributions  to  these  displacements 
arising  from  (a)  thermal  expansion,  (b)  axial  forces,  (c)  shear 
forces,  and  (d)  bending  moments.  The  moments  are  related  to 
the    axial    forces    by   the    moment    equilibrium    equation    which    gives 


\    2  Tx  (x)    +       2       3  T2  (x)    =  M1   +  M2   *  Mz 


(7) 


where  T,  and  T2  are  the  axial  forces  on  layers  1  and  2.  The  Tj 
axial  forces  are  related  through  equilibrium  to  the  shear  stresses 
between    layers,    that    is 


T^x)   =  fxtiZ)  dl 


(8) 


where  r^x)  is  the  shear  stress  between  the  ith  and  (i  +  1)st 
layers.  The  M,  moments  are  related  to  the  radii  of  curvature,  ps 
through 


Md(x)    =  -  — %- 
p  (x) 


i   =  1,2,3 


(9) 


and    the    D,    flexural    rigidity    coefficients    for    each    layer    are    given 
by 


*>i  = 


E±hl 


12{l-\\) 


(10) 


where  E„  hj,  and  v,  are  the  Young's  modulus,  thickness,  and 
Poisson's  ratio  of  the  ith  layer  respectively.  Equation  (10)  assumes 
there  is  no  variation  of  lateral  displacement  with  respect  to  y, 
and    therefore    the    radii    of   curvature    of    layers    are    all    equal. 


Using    equilibrium   equation    (7)    and    moment-curvature   equation    (9) 
yields 


pU) 


=   -Tx(x) 


-  r,  (x) 


A. 


12 


(ID 


where 


D  =  dx  +  D2  *  D3 


(12) 


is    the    flexural    rigidity    of   the    assembly,    and    where 


•i,i*l 


(t>i  +  hi.iY 


2D 


(13) 


is     an     additional     compliance     due    to     bowing.         The     individual 
contributions    to    Eqn's    (5)    and    (6)    are   the 


•  the    thermal    contribution    is: 

^(AD^  (14) 

•  the    axial    force    contribution    is: 

-kA'TiWdt  (15) 

J  o 

where  Xy   is  the  axial   compliance  of  the   i,h   layer,   given   by 


•    the    shear    stress    contribution    is:    -/c,t,(x),    where    k-{    is    the 
interfacial    compliance    given    by 

K,    =    — i-  (17) 

1         3G, 


•    the    rotational    contribution    is: 

-hifx    dl 


i 


2    Jo     p(£) 


(18) 


Invoking    the    condition    that    the    shear    stresses    be    anti-symmetric 
with    respect   to   the   origin    of   coordinates   leads   to   the   expression 

xt  (x)   =  Cjsinh  kx  i=if2  (19) 


where  k  is  obtained  by  solving  a  characteristic  equation  defined 
in  Eq.(20)  below.  Substitution  of  the  above  expressions  into 
compatibility    conditions    (5)    and    (6)    yields   the    quartic    equation 

K12K33^4-  <Kia*23+K23Au>  k'    +    *    =    0  (20) 

where  the  kv's  and  Ay's  are  interfacial  and  axial  compliances, 
respectively,    of   the   tri-material    assembly,    given    by 

K12    =   K,    +    K2  (21) 


8 


*23  *2     +    K2 


(22) 


o 

'23 


X2  +  A3  +  A23 


(23) 


A12 


=  *i   +  *2   +  *i: 


(24) 


6  =  xx>.2  +  x2x2  +  x3x,  +  a3a12  +  x2(/xt;  +  yx^) 

Letting    to    =    k2,    Eq.(20)    becomes 


K12K33G)2     ~     (K12*23     +    *23X?2)W     +    6     =    0 


whose    solution    is 


GO 


1,2 


-b±Jb2  -4ac 


2a 


(25) 


(26) 


(27) 


where 


a  -  k12k?3 


b  =  -(k12a23  +  k23*i2) 


The    k,    (i  =  1,2)    roots    are   then    obtained    as 


(28) 


(29) 


The     expressions      presented      above,      valid      for      any     tri-material 
assembly,    reduce    to    the    simpler    form 


T(x)    =  k  ,AgA/   sinhicx 
Acosh  1 


(30) 


where 


k  = 


^ 


0 
■2? 


K    =   Kx    + 


2^    +    K3 


(31) 


when  the  Youngs  modulus  of  the  mid  layer,  or  its  thickness  is 
very  small,  as  in  the  case  of  a  thin  layer  of  solder  used  as  an 
adhesive.    In    accordance   with    Eq. (1 9),    the    maximum    shear    stress 

occurs    at    x    =    f,    and    can    be    calculated    from 
where 


■  *^M^  <"> 


.  i        cosh  k{ 1-a)  ,--» 

X~*  "  1  ~         cosh  kl  {      } 


Application    of    Suhir    Model    to    Unleaded    Device 

Here,  the  closed  form  equations  from  Suhir's  model  are  used  to 
determine  the  shear  stresses  in  an  uleaded  device.  The 
application  is  to  unleaded  device  with  dimensions  and  material 
properties    as    follows: 

Material    #1;    Ceramic    chip: 

a,     =    3.6x106    in/in/F,    v    =    0.3,    E,     =    40x106    psi, 
G,     =     15.4x106    psi,    h,     =    .056    in. 

Material    #2;    Solder: 

a,    =    25.X10'6    in/in/F,    v    =    0.3,    E2    =    1.86x106    psi, 
G2    =    .66x106    psi,    h2    =    .002    in.  ' 

Material    #3;    (board): 

a3    =    16.x10'6    in/in/F,    v    =    0.4,    E3    =    2.0x106    psi 
G3    =    .77x106    psi,    h3    =    .056    in.  ' 

The    above    parameters    lead    to   the   following    compliances: 

Chip:  D,    =    643    lb. in.,         X,    =    .3125x10"6    in. /lb., 

k,    =    1.212x10'"    in. /lb. 

Solder:    D2    =    .0231    lb. in.,    X2    =    64.52x10"6    in. /lb., 
k2    =    2.525x10"9    in. /lb. 

Board:       D3    =    32.2    lb. in.,      A3    =    6.25x10'6    in. /lb., 
k3    =    24.24x10-"    in. /lb. 

The    other    parameter   values    necessary    for    a   solution    are: 

A12    =    1. 378x1  0'6    in. /lb.  A23    =     1. 378x1  0-8    in. /lb. 
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/c12    =    3.737x109    in.  /lb.  k2  3    =    26.77x109    in.  /lb. 

A12'    =    64.83x10"6    in. /lb.  A23'    =    70.77x106    in. /lb. 

A12°    =    66.21x106    in. /lb.  A23°    =    72.15x106    in. /lb. 

Substitution    of   these    parameters    into    Eq.(25)    yields 

b    =    .7897x109    (in. /lb.)2 

which    upon    substitution    into    Eq.(26)    yields    k    =    19.85    in.1. 
Finally,     the     maximum     shear     stress     at     the     solder     interface     is 
obtained    from    Eq.(32)    as 

TmaY    =    4360    psi 

max  " 

The  simpler  approximate  expression  for  k  from  Eq.(31)  yields  k 
=     19.61     in.'1,    which,    in    turn,    gives    rmax    =    4310    psi. 

Finite    Element    Model 

An  alternate  model  of  the  tri-assembly  problem  using  a  finite 
element  method  was  undertaken.  Continuity  along  the  interfaces 
of  adjacent  layers  is  maintained  by  selection  of  the  six  degree 
of  freedom  element  shown  in  Figure  8.  The  six  degrees  of 
freedom    are    as    follows: 

d,  is  the  axial    displacement    at    the    lower    left    corner 

d2  is  the  axial    displacement    at    the    upper    right    corner 

d3  is  the  lateral    displacement    at    the    left    end 

dA  is  the  axial    displacement    at   the    lower    right    corner 

<35  is  the  axial    displacement    at   the    upper    right    corner 

66  is  the  lateral    displacement    at    the    right    end 

The  axial  displacement  field  u(x,y)  is  taken  as  linear  interpolation 
in    both    the    axial,    x,    and    transverse,    y,    directions.    That    is, 

u{x,y)  =  H1(y)  (Nx  (x)  6X  +N2  (x)  62 )  +H2(y)  (Nx  (x)  64  +N2  (x)  65)     (34) 
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and     the     lateral     (i.e. .transverse)     displacement     field     is     a     linear 
function    of    x,    that    is, 


.   ,  «.!     .    \*   W  «    W    , 


v(x)    =  ^(x)  63   +  N2(x)  66  (35) 


where    the    linear    interpolation    functions    are: 


Hx{y)    =  ±(h-y)  (36) 

n 


H2(y)    =   Z  (37) 


Njx)    =   j  (1-x)  (38) 


NAx)    =   4  (39) 


"2 


1 


The    strain-displacement    relations    are: 


.,--i-^<«,-»x>*^- <».-»,>         «°> 


^f'I  =  ^<5^>^<^>^  ««> 


On    the    element    level,    we    have   the    bending    effects 

Mb  "  Fb  <*2) 

where 

S£  =  <61    62    64    65>  (43) 

and 

KB  =  f1fhBEBTdydx  (44) 

Jo  Jo 
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and 


B 


T    _ 


d^  dNx  dN2  dN2 

<"&T/fl    "a^"2    ~teHi    ~teH^ 


(45) 


and 


FB  =   f    [   BEBTdydx 
Jo  Jo 


(46) 


which    gives 

r  2 

1 

-2 

-1 

B         6  1 

1 
-2 

2 
-1 

-1 
2 

-2 
1 

L  1 

-2 

1 

2 

Similarly    for    the 

shear    effects 

K 

?6s 

=   0 

(47) 


(48) 


where 


65  =  <br    62    63    64    65    66  > 


(49) 


and 


Ks  =    f    fhB'GBITdydx 
Jo  Jo 


(50) 


where 


,T  dH,  dH7       dN,  dH,  dH7       dN7 

B       —   <  N 


N- 


dy      "1  dy        dx        2  By        2  dy        dx 


>  (51) 


which    gives 
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Ks  = 


-a 

G 

2 


-a 


G 
2 


-a       — 


a     -  —     -a 


-  a 


-iZ        b      J2     -if     -jb 


■a      — 


-—     -a 


—  -h      -  _ 


G 
2 


G 
2 


G 
2 
G 
2 


-a    -  — 


a      — 


G 
2 

G 
2 


(52) 


The  bending  and  shear  element  stiffness  matrices  are  combined 
to  form  the  element  stiffness  matrix,  which  in  turn  are  again 
combined  to  form  the  system  stiffness  matrix.  In  accordance  with 
standard  FEM  techniques,  the  system  matrix  equations  are  solved 
for  the  system  displacement  degrees  of  freedom.  The  element 
strains  are  then  calculated  from  Eq.(40)  and  (41)  above.  Finally, 
the  average  bending  stresses  at  the  bottom  and  top  of  an 
element    are    obtained    from    the    stress-strain    relations 


(°x>l 


.      »4-*l 


(53) 


and 


oJ,  = 


'  x>  2 


_  s5-^ 


(54) 


Similarly,    the   average    element   shear   stresses   are   calculated    from 
the    strress-strain    relations 


(Oi  = 


62_61  66_53 


(55) 


and 


^xyh    = 


65-64        66-6. 


(56) 
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The  finite  element  mesh,  shown  in  Figure  9,  consisting  of  104 
elements  and  149  degrees  of  freedom  produced  an  interface 
maximum  shear  stress  of  3650  psi.  This  FEM  calculated  stress  is 
at  the  midpoint  of  the  right  most  element.  This  stress,  obtained 
from  an  FEM  analysis  is  at  the  midpoint  of  the  interface  at  the 
right  most  element.  Quadratic  extrapolation  to  approximate  the 
hyperbolic  sine  expression  for  the  r  shear  stress  brings  the  FEM 
predicted  value  to  about  3800  psi.  Thus  the  predicted  values  of 
the  maximum  shear  stress  from  the  Suhir  model  and  the  FEM 
model    differ    by    less   than    15    percent. 
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Conclusions 

This  investigation  has  shown  that  under  the  same  thermal 
environment,  the  stiff  unleaded  solder  joint  has  significantly  higher 
thermoelastic  stresses  than  the  flexible  leaded  solder  joint.  Since 
the  unleaded  joint  has  a  stronger  suppression  of  displacement 
than  the  leaded  joint,  and  thermal  stresses  result  from 
displacement  constraint,  the  result  is  not  an  unexpected  one.  In 
the  particular  application  presented  in  the  previous  sections,  the 
stresses  in  the  unleaded  connection  were  twenty  times  the 
stresses  in  the  leaded  connection.  This  large  degree  of  difference 
between    the    stresses    in    the   two    cases    is    surprising. 

As  previously  noted,  the  comparison  is  for  the  case  where  both 
connections  are  subjected  to  the  same  thermal  environment.  This 
assumption  does  not  take  into  account  the  fact  that  the  thermal 
field  is  very  likely  to  be  different  for  the  two  types  of  joints. 
One  might  expect  that  the  leaded  joint  with  the  'separated' 
components  may  have  larger  thermal  gradients  than  the  unleaded 
joint  where  the  components  are  in  close  proximity  to  each  other 
and  therefore  may  have  a  more  uniform  thermal  field.  In  fact,  the 
solution  to  the  stress  problem  really  requires  its  coupling  to  the 
thermal  (energy  balance)  problem.  This  later  coupled  (and 
transient)  problem  is  presently  under  investigation  as  a  IMPS 
thesis   topic   under  the   direction   of   D.    Salinas   and   Y.W.    Kwon    (2). 
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Figure    2.    Cross    section    of    lead    wire 
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Figure    5.    Deformation    of    beam    BC 
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Figure    8.    Six    degree    of   freedom    element 
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Figure    9.    Finite    element    mesh 
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